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Randall-Sundrum graviton spin determination using azimuthal angular dependence 
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Quantum interference of helicity amplitudes provides a powerful tool for measuring the spins of new particles. 
By looking at the azimuthal angular dependence of the differential cross-section in the production followed by 
decay of a new particle species one can determine its spin by looking at the various cosine modes. The heavy 
spin-2 Kaluza-Klein (KK) graviton provides a unique signature with a cos (4<^) mode. We study the feasibility 
of this approach to measuring the spin of the KK graviton in the Randall-Sundrum Model at the LHC. 

PACS numbers: 



I. INTRODUCTION 

The Large Hadron Collider (LHC) at CERN is expected to 
produce a wealth of discoveries by probing the TeV scale for 
the first time. Apart from finally accessing the electroweak 
symmetry breaking scale and thus potentially discovering the 
elusive Higgs boson, we expect to see new physics that re- 
solves the hierarchy/naturalness problem |[ll 0. B 0] and per- 
haps provides an insight into the nature of dark matter. One 
exciting possible solution to the hierarchy problem is the exis- 
tence of warped extra dimensions |5, 6] which allows for TeV 
scale gravitational interactions. There are many variations of 
the basic theory [7] but one common feature that they share is 
the existence of heavy Kaluza-Klein (KK) gravitons. 

For the purpose of this paper we will consider a Randall- 
Sundrum model with 3h-1 dimensional spacetime with one ad- 
ditional warped extra dimension (RSI). The Standard Model 
fields are confined to a 3h-1 dimensional TeV brane and the 
graviton propagates freely in the 4h- 1 dimensional bulk. Quan- 
tization of the graviton wave function in the extra dimension, 
with boundaries between the TeV brane and a Planck brane, 
leads to various modes which appear as heavy spin-2 fields in 
the 3h-1 dimensional effective theory on the TeV brane. 

One of the challenges at the LHC will be to determine the 
spins of newly discovered particles in order to distinguish var- 
ious theoretical models. The KK graviton provides a unique 
signature of gravitational physics at the TeV scale by virtue of 
its spin-2 nature. Thus, it becomes crucial to have techniques 
to identify its spin. So far, the technique proposed to measure 
KK graviton spin at the LHC relies on resonant graviton pro- 
duction followed by decay into a lepton pair [8, 9, lO. illlll^ . 
By looking at the polar angular dependence of the leptons rel- 
ative to the beam axis, one sees a quartic behaviour of the 
differential cross-section. 



da 
decs 6* 



Acos^e + Bcos^e + G 



(1) 



Recently, a new technique for measuring spin has been 
proposed. One can look at quantum interference of helicity 
states in the azimuthal angular dependence of particle decays 
to study their spin in a model independent way lITsilTill . The 
goal of this paper is to apply this technique to study the KK 



graviton spin and look at its feasibility at the LHC. 



II. MODEL PARAMETERS 

The interaction between the massive KK gravitons and the 
Standard Model fields in the 4-d effective theory is given by 
the Lagrangian lfT6l[T7ll . 



(2) 



Here, G^^^'^'^ represents the nth KK graviton mode. 7^;^ is the 
stress-energy tensor of the Standard Model Lagrangian given 

by. 



A is the coupling given by. 
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where k is of the order of the Planck scale, Tc is the compact- 
ification radius of the extra dimension and Mpi = Mpi/y/Sn 
is the reduced Planck scale. Note the absence of KK-parity 
which allows the heavy graviton modes to decay into purely 
Standard Model particles. 

The mass of the ?ith KK-graviton is given by. 



XnA- 



(5) 



where a;„ are the nth zeros of the Ji Bessel function. While 
studying the properties of the n = 1 KK graviton we can thus 
regard this theory as being dependent on only two parameters 
A and k or equivalently the dimensionless coupling c = -ry— 
and nil, the mass of the KK graviton of interest. 

Naturalness constraints require A less than about 10 TeV. In 
order for an effective field theory description of gravity to be 
valid we require that the 5-d curvature bound, |i?5| < A/|, is 
satisfied, where A/5 is the 5-d Planck scale. By looking at the 
various theoretical and experimental constraints on the model 
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FIG. 1: Experimental and theoretical constraints on the KK gravi- 
ton parameters in the c — mi plane. Red curves show experimental 
constraints and blue curves show theoretical constraints. The green 
shaded region shows the allowed parameter space. 



parameters (Figure [Hi we expect c to lie roughly between 
0.01 (weakly coupled) and 0.1 (strongly coupled). We con- 
sider TOi in the range of 750 GeV - 2 TeV. The decay width of 
the graviton to Standard Model particles can be evaluated by 
using the expressions given in LIQ, JLT, J^8J. In the limit that 
decay particle masses can be neglected the decay width of the 
graviton is given by 



am, 



(6) 



where a is a constant depending upon the number of open 
decay channels. If one assumes decay to only Standard Model 
particles the ratio Fi : nii is found to be 1.37% for c = 0.1 
(Assuming a Higgs mass of 120 GeV and decay into Standard 
Model particles only). This value is in disagreement with the 
value 1.43% cited in the literature I,l9i1 . 



III. USING AZIMUTHAL ANGULAR DEPENDENCE TO 
MEASURE SPIN 

To determine the spin of a particle X, we consider the pro- 
duction process A + B X + Y where X further decays to 
M + N. Here, A and B refer to beam particles or partons, X 
is the parent particle whose spin we wish to measure. M and 
N are the daughter particles that X decays into. 

This gives us two planes to consider, namely the production 
plane (defined by the beam direction and the parent momen- 
tum) and the decay plane (defined by the parent momentum 
and either daughter) (Figure |2]l. 

Now consider the daughter M with momentum pM- The 
angle it makes with the parent momentum px is defined to 
be 9. Projecting out the component of pM parallel to px and 
looking at the angle between the residual vector and the pro- 
duction plane we define an angle 0. Thus, (j) describes az- 
imuthal rotations of the vector pM in the x — y plane with px 
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FIG. 2: Production and decay planes of the process A + B ^ X + Y 
AI + N. The angle (j> is defined as the azimuthal angle between 
Px and pAf or equivalently the angle between the production and 
decay planes. 



taken to be the z-axis. From the figure it is clear that, equiv- 
alently (f) can be defined as the angle between the production 
plane and the decay plane. More explicitly, we define the two 
vectors. 



and. 



Then, 



Pa X Px 



Pdecay = X pM 



Pprod X Pdecay 



(7) 



(8) 



(9) 



Here p denotes the normalized vectors. 

In the limit of the narrow width approximation, the ampli- 
tude can be split into Mprod and M decay 



M 
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{X,Y\Tprod\A,B) 



M decay 

{cl,)^{M,N,cf>\%ecay\X) 



(10) 



(11) 



where we have explicitly shown the (j) dependence of the the 
final state and decay amplitude. We also have. 



MdecayW = {M , N = 0)\e+''-''^%e 



\X) 



(12) 



where generates rotations about the px direction. We can 
now think of the rotation operator as acting on the interaction 
T-matrix plus ket, rather than on the bra. Assuming, Tdecay 
is rotationally invariant, we only need to consider rotations of 
the particle X about its own momentum axis. In this case. 



Jz — J ■ p — {s + fx p) ■ p ~ s ■ p 



(13) 



Thus, rotations about the momentum axis of a given helicity 
state, h for X only produce a phase e+*'"^. So, 



Mdecayi(l)) = M decay = 0) 



(14) 



Thus, allowing for production over all possible heUcity states 
of X we must sum coherently over all possible amplitudes, 
and so, the differential cross section takes the form 



da 



Y,^Prode+''"^Mdecay{(l> ^ 0) 



(15) 
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Here, h runs from —s to +s where X has spin s. From this it 
is clear that, if we look at the differential distribution da/dcj), 
interference between various helicity states is responsible for 
a non-trivial (p dependence, 

^ ^ Ao+Ai cos (0) + A2 cos (20) + .... + cos (2s0). 
d(p 

(16) 

Note the absence of sin {mp) modes, which would be present 
in the case of CP violating processes. 

The Standard Model has no particles with spin greater than 
1 and so the largest mode from the Standard Model would 
only be cos (20), corresponding to X being a gauge boson. 
We can now see the unique signature that the KK graviton 
will produce, namely a cos (40) mode. 

Also, we note that this result is valid in any reference frame 
but the size of the coefficients Ai will be different in different 
reference frames. To maximize this unique signature for the 
KK graviton, we need to choose a reference frame in which 
A4/A0 has a large value. 




Lab Frame/CM Frame Zero-Rapidity Frame 



FIG. 3: Boost from center of mass or laboratory frame to the Zero 
Rapidity frame 

large. It was found that in the center of mass frame of the par- 
tonic processes, S4 was larger than in the lab frame. However, 
transforming from the lab-frame to the center of mass frame 
would have an error dependent on the error of the jet recon- 
struction. To avoid this error and still make an improvement 
in the signal, we studied ^4 in the zero rapidity frame of the 
graviton, i.e. the frame where the graviton is purely transverse 
to the beam axis (Figure [3]). The reason for this is the boost 
factor can be calculated from just the graviton momentum in 
the lab frame which is well reconstructed from the leptons. 



IV. SIGNAL AND BACKGROUND 

We assume that the mass of the graviton will be well mea- 
sured using resonant graviton production through the process 
pp^G^l+l- |l8t|9|]. 

The process we are considering is ^ G H- jet followed by 
G ^ where I are muons or electrons. The dominant par- 
ton level subprocess comes from gg Gg with subdominant 
9(9) 9 ^ G q{q) and the crossed channel q q ^ G g. Here G 
represents the graviton, g represents gluons and q represents 
the various quarks. 

The Standard Model background comes from the subdom- 
inant channels with G replaced by an off-shell Z, 7. This is 
the exact analog of Drell-Yan background in resonant graviton 
production. Cutting on the invariant mass of the lepton pair in 
a mass window around the graviton mass gets rid of most of 
the background. The Standard Model background consists of 
spin-1 states and can not give any contribution to A4. At most 
it can affect the value of and dilute the value of ^4/^0- 



V. CALCULATING THE DIFFERENTIAL 
CROSS-SECTION 

A. Zero Rapidity Frame 

The dilepton + jet events that we are looking for are fully 
reconstructible at the LHC. The key reason for this is that we 
have a signature with no missing energy-momentum which in 
turn is a direct consequence of the absence of KK parity. The 
graviton 4-momenta should have minor errors compared to the 
jet reconstruction since it is reconstructed from the di-lepton 
4-momenta. As previously mentioned the size of the non-zero 
coefficients Ai are frame dependent and so we must choose a 
frame in which the normalized coefficient S4 = 1^4/^0! is 



B. Cuts 

The first set of cuts used included a pseudo-rapidity (I77I < 
2.5) cut and pT > 20 GeV cut for the jet. The second set of 
cuts was a mass-window cut on the invariant mass distribution 
of the lepton pair. This gets rid of a large portion of the Stan- 
dard Model background. The size of the window was deter- 
mined by detector resolution at ATLAS lf22ll23ll for an e+e^ 
pair The third set of cuts involved rapidity cuts {\ri\ < 2.6) 
on each of the leptons with a requirement that pT > 10 GeV 
for either one of the leptons and px > 20 GeV for the other. 

An isolation cut, Ar = {Arj)'^ + A(j)^ > 0.7, was imposed 
between the lepton and the jet. However, the third set of cuts 
affects the angular distribution of the leptons and can create 
false cosine modes in the differential distribution. 

To solve this problem one imposes 'rotationally invariant 
cuts', first introduced in (T^. Thus, it is not sufficient for the 
observed lepton to simply pass these cuts, the leptons are ro- 
tated around the graviton momentum axis in small increments 
and at each step it is checked that the lepton passes the cuts. 
The added complication is that the rotations must be made 
in the zero-rapidity frame to preserve rotational invariance in 
that frame (FigureHj). So, the procedure is as follows: 

1 . First, reconstruct the event completely using the dilep- 
ton and jet signals. 

2. Calculate the boost factor to take us from the lab frame 
to the zero-rapidity frame of the graviton. 

3. Boost all momenta to the zero-rapidity frame. Rotate 
the leptons about the graviton momentum direction by 
a small angle, say 1°. 

4. Reboost the new lepton and jet momenta to the lab 
frame. Check if they pass the cuts, if they don't throw 
out the event. 
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Incremental Rotations 



Dilepton + Jet signal 




FIG. 4: The leptons are rotated about the graviton momentum axis 
in the zero rapidity frame. The dilepton + jet momenta must be re- 
boosted to the lab frame at each step to make sure that they pass the 
cuts. 
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1 TeV and c = 0.05. 
A strong cos (2<j!>) mode can be seen but there is also a cos (4<^) com- 
ponent. The theoretical curve (produced from simulations) is shown 
in green. The red dots indicate the binned values, with error bars 
corresponding to Gaussian errors for a luminosity of 500 fb^^ 



5. If they do pass the cuts go back to step 3. 

6. Repeat this procedure until we have made a full 360° 
rotation of the lepton momenta about the graviton mo- 
mentum axis in the zero rapidity frame. 

This procedure ensures that the cuts do not affect the az- 
imuthal angular distribution in the zero-rapidity frame. 



C. Simulations 

We used HELAS |24] with spin-2 particles lEsIl to calcu- 
late the helicity amplitudes for the graviton scattering process. 
LHApdf ll26ll was used to fold in the parton distribution func- 
tions for the protons. We used the pdf set CTEQ6L Il27[l . An 
adaptive Monte-Carlo package, BASES [281, was used to per- 
form the integration over phase space and produce the differ- 
ential cross-section da/d(t). 



where i runs over 0, 1, 2, 2n—l. The integration accounts 
for the binning process and the 2?! coefficients Aq, An-i, 
Bi, Bn correspond to the strengths of the various cosine 
and sine modes that can be resolved from each other. 

Thus, we have a simple linear relationship between the 2n 
binned values of da/d4) (xi) and the 2n binned-Fourier coef- 
ficients (j/j) of the form Xi — PijUj- Here, pij are either of the 

27vi 2Tvi 

form / cos {j(j))d(l> or / sin {j(j))d(l>. 

Now, we can simply invert this matrix for a given value of 
n to recover the amplitudes of the various harmonics. For the 
da/ dip distribution for the graviton we expect to see only the 
coefficients Aq, A4 to be non-zero. Also, since the beams 
are identical, we expect to see only the even cosine modes. 
The odd cosine modes drop out since they flip sign when the 
beams are switched (0 ^ vr — (p). 



VI. DETERMINING THE COEFFICIENTS OF THE 
VARIOUS COSINE MODES 

Once we have the binned distribution da/dcj) with 2n bins 
(For the purposes of calculation in this paper we used 50 bins), 
we try to fit coefficients of the form 



1 



Binsize / ^^(i-i 



da 



27r/2n Jliitill 



Aj cos (j0) 



+ Y,Bj sin {j4>) 



(17) 



VII. RESULTS AND DISCUSSION 



Simulations were done for the process pp 



J at 7 

TeV beam energy using a dilepton invariant mass window cut 
around the graviton mass. Figure|5]shows the da/dcp distribu- 
tion for a 1 TeV graviton with c — 0.05. Figure |6] shows the 
normalized fitted coefficients. The normalized cosine coeffi- 
cients (Si) are shown in the first 25 bins, with the zero mode 
suppressed. The next 25 bins show the sine modes. The size 
of the 54 coefficient is 3.14%. Note the absence of odd cosine 
modes, this arises from the fact that we are using identical 
beams. 

To look at the dependence of the signal on graviton mass, 
simulations were performed for c = 0.1 and mi = 750 GeV, 
1 TeV, 1.5 TeV and 2 TeV. The results are summarized in Ta- 
ble U The total cross-section decreases rapidly with graviton 
mass as expected. The background is negligible and, as we 
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FIG. 6: Fitted cosine coefficients of the binned differential cross- 
section shown in Figure |5] corresponding to 50 bins. The first 25 
modes label the normalized cosine modes, the next 25 show the sine 
modes. The large 0-mode which would be 100% is not shown). See 
text for how the error bars in this plot are calculated using error bars 
from the binned differential cross-section. 



mi (TeV) 


Am (GeV) 


CTtotal (fb) 


O-bgd (fb) 


S2 


S4 


0.75 


24.4 


871.7 


0.39 


20.00% 


3.50% 


1.0 


30.7 


229.8 


0.15 


20.48% 


3.16% 


1.5 


42.8 


28.7 


0.03 


20.70% 


1.52% 


2.0 


55.0 


5.52 


0.01 


20.08% 


0.80% 



TABLE I: Signal strength ^4 = | /ylo | as a function of the mass of 
the graviton. c = 0. 1 for all entries. ^2 is shown for comparison. The 
mass window (based on the ATLAS detector resolution for e^e~ 
invariant mass I23.l23ll ) cuts out most of the background. 



will see in the next paragraph, has little effect even if the cou- 
pling c is reduced. The main concern is therefore the decrease 
in 5*4 and the low cross-section at large values of mi. 

The results for a 1 TeV graviton at different values of the 
coupling c are shown in Table [III In the absence of cuts the 
graviton cross-section is expected to approximately scale like 
c^. The Standard Model background level is 0.15 fb which 
is ~ 5% at c = 0.01. The value of S4 is expected to be di- 
luted slightly by the background because of a corresponding 
5% increase in Aq. As c is increased, the background as a 
percentage of the cross-section decreases and S4 is restored to 
its maximum strength. 



VIII. ERROR ANALYSIS 



As we have seen, the effect of background is small and 
does not contribute to A4. Its only effect is to dilute the nor- 
malized coefficient S4. Thus, the experimental error will be 
determined by event statistics. We assumed Gaussian errors 



(Axj — Xj ^j^"' ) in the jth bin assuming Nj = ^'^^ — : 

events in each bin for integrated luminosities £ of 10, 100 and 
500 fb~^. Since, the coefficients Ai are determined from the 
binned values xj through a simple linear relationship (via the 
matrix = PTj^^ '■^^'^ straightforward to work out the 
errors in the normalized coefficients (ASi). 



AS, = 



\ 



j 



Ax? 



An 



A, 



(18) 



The first term in the paranthesis arises from the simple lin- 
ear relationship between A.^ and xj . The second term comes 
from the eiTor associated with the normalization factor Aq. 
The relative errors (^^) for various integrated luminosities at 
different points in the parameter space of the model are given 
in Table |lll] A value > 1 for the relative eiTor indicates that 
statistics would be poor and give no reason to doubt 54 being 
consistent with 0. A value of 0.20 or less indicates at least 
a 5(7 effect indicating high likelihood of confirmation of the 
spin-2 nature of the KK graviton. 



mi (TeV) 


c 


lOfb"^ 


100 fb"^ 


500 fb-^ 


0.75 


0.1 


0.43 


0.14 


0.06 


1.0 


0.01 


8.03 


2.54 


1.14 


1.0 


0.02 


3.97 


1.26 


0.56 


1.0 


0.05 


1.65 


0.52 


0.23 


1.0 


0.1 


0.93 


0.29 


0.13 


1.5 


0.1 


5.42 


1.71 


0.77 


2.0 


0.1 


23.52 


7.44 


3.32 



TABLE III: Statistical Error AS'4/5'4 for different integrated lumi- 
nosities for the process pp e^e~j. AS4/S4 < 0.5(0.71) cor- 
responds to a 2(7 confirmation of the graviton spin, and AS^/S^ < 
0.2(0.28) corresponds to a 5a confirmation. The values in brack- 
ets denote the 2a and 5a confidence levels if one includes n'^ fj,~j 
production channels as well. 



c 


atotai (fb) 


S2 


S4 


0.01 


3.27 


18.62% 


3.05% 


0.02 


12.51 


20.02% 


3.15% 


0.05 


72.75 


20.42% 


3.14% 


0.1 


229.8 


20.48% 


3.16% 



TABLE II: Signal strength S4 = \A4/Ao \ as a function of the cou- 
pling c. All entries are for mi = I TeV. S2 is shown for comparison. 
The Standard Model background cross-section is 0.15 fb. 



Alternatively, if one requires only a 95% confidence level 
(2a-) effect then a value of 0.5 or less for AS'4/S'4 should 
suffice. If we additionally assume information from fi^fi^j 
statistics in addition to the e+e~j channel (assuming that de- 
tector resolution for the invariant mass is the same for both 
lepton species), then we can see a factor 2 improvement in the 
statistics. This would in turn result in a factor \/2 drop in the 
error. Thus, in this case the parameter space in Table Hill for 
which AS'4/S'4 < 0.71 would correspond to potential for a 2a 
confirmation of the graviton spin, and AS'4/S'4 < 0.28 would 
correspond to a 5cr confirmation. 
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IX. COMPARISON WITH RESONANT GRAVITON 
PRODUCTION METHOD AND DISTINCTION FROM 
SPIN-0 

Osland et al. ifTll [T2I1 consider the resonant graviton pro- 
duction process pp G l^l^ to measure the spin of 
the graviton using the quartic angular dependence of the polar 
angle of the lepton. This results in a center-edge asymme- 
try {Ace) in the differential distribution da /d cos 9. Their 
results indicate {2a) identification of the graviton spin for 
c — 0.01 and 10^^ fb of luminosity for masses upto 1.1 TeV. 
For c = 0.1 they claim identification upto masses of 2.4 TeV. 

The azimuthal angular dependence method that we con- 
sidered has inherently lower statistics compared to resonant 
graviton production because of the extra recoiling jet. Our 
method suffers from lower statistics, but given higher lumi- 
nosities, it can still provide an independent confirmation of 
the KK graviton spin for a large region of the expected pa- 
rameter space of the KK graviton. 

The center-edge asymmetry method can distinguish a spin- 
1 particle {Z') from a KK graviton more readily than it can 
distinguish it from a spin-0 particle. 

Our method proves complementary, since the KK graviton 
also produces a large cos(20) mode (^2 ^ 20%) and can thus 
easily be distinguished from a scalar which would not produce 
any non-zero modes. The results for AS'2/S'2 are shown in 
TableHV] 



mi (TeV) 


c 


10 fb"^ 


100 fb"^ 


500 fb"^ 


0.75 


0.1 


0.07 


0.02 


0.01 


1.0 


0.01 


1.30 


0.41 


0.18 


1.0 


0.02 


0.62 


0.19 


0.09 


1.0 


0.05 


0.25 


0.08 


0.04 


1.0 


0.1 


0.14 


0.04 


0.02 


1.5 


0.1 


0.39 


0.12 


0.06 


2.0 


0.1 


0.93 


0.29 


0.13 



TABLE IV: Statistical Error AS2/S2 for different integrated lumi- 
nosities for the process pp — > e'^e~j. A.S2/S2 < 0.5(0.71) cor- 
responds to a 2(7 distinction from a spin-0 particle, and A5'2 / S2 < 
0.2(0.28) corresponds to a 5cr distinction. The values in brackets 
denote the 2a and 5a confidence levels if one includes ji^ ^~ j pro- 
duction channels as well. 

Assuming, as before inclusion of fJ.^fJ,^ statistics 
AS2/S2 < 0.71 corresponds to a 2a distinction from a spin-0 
particle, and AS'2/S'2 < 0.28 corresponds to a 5a distinction. 
In regions of the parameter space of mi, where the Standard 
Model background is comparable to the cross-section of inter- 
est (TablelH, the confidence levels are altered slightly because 



the off-shell 7 and Z, being spin-1, contribute to the A2 coef- 
ficient. 

From the table, we can see that even with 10 fb^^ of lumi- 
nosity, the spin-0 hypothesis can be ruled out for a large por- 
tion of the allowed parameter space. Thus, our method proves 
complementary to the approach by Osland et al. by ruling 
out spin-0 more easily than spin-1. In both methods the dis- 
tinction from spin-0 can be made from comparable integrated 
luminosities. 



X. SUMMARY AND CONCLUSION 

We studied the process pp ^ G jet ^ jet and looked 
at the differential distribution da/d(p. The distribution was 
found to have a cos (4(/)) mode, characteristic of a spin-2 par- 
ticle, with strength parametrized by S4. The parameter S4 
was ^3% for values of toi below a TeV. As we go to higher 
graviton masses the signal drops off, but what is of more con- 
cern is the drop in cross-section with large mi or low values 
of c. Both these scenarios are unlikely to occur in conjunction 
because of naturalness constraints (See Figure[T]i. 

In conclusion, observing higher cosine modes (> 2) in 
the differential distribution would be a clear signal of beyond 
Standard Model physics. Observing the cos(4(/)) mode at the 
LHC would be a strong indicator of gravitational physics at 
the TeV scale. If the coupling is strong enough ^ 0.05 or 
greater and the mass is sufficiently low ^ 1 TeV or less, we 
expect to have a clear signal of the spin-2 nature of the KK 
graviton at the LHC. 

For regions of parameter space with larger masses or lower 
couplings, the azimuthal angular dependence of the cross- 
section is still useful in ruling out the spin-0 hypothesis and 
this can be done for fairly low luminosities ^10 fb^^ as well. 

This method provides an important complementary and in- 
dependent approach to measuring the spin of the KK graviton, 
as compared to the method of using polar angular dependence 
from resonant KK-graviton production. 
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